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ABSTRACT: 

An analysis of a 3- inch gun launched finned motor case is performed in an 
attempt to determine the cause of failure. Two specific features of the problem 
are identified as the probable reasons for the buckling failure of the motor case. 
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INTRODUCTION 



The Naval Weapons Center, at China Lake, California, is presently 
in the process of developing a series of gun launched guided projectiles. 
Recent experimental firings of a 3-inch projectile showed the structural 
design of the motor case to be inadequate since the case buckled severely 
as the projectile passed through the gun barrel. The authors were asked 
to perform an analysis of the motor case to determine the cause of 
failure. In addition, the authors were asked to recommend a structural 
design of a 5 inch projectile that minimizes the motor case weight and 
to propose a design technique which can be used by engineers as a 
design tool in future gun launched guided projectile projects. This 
report concentrates on the analysis of the 3-inch motor case and the 
design technique. 

The following four pages show l) a drawing of the projectile; 

2) and 3) high speed photographs of two firings; and 4) a picture of 
a buckled projectile in flight. 
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DESCRIPTION OF THE THREE-INCH PROJECTILE 



The casing is 4130 steel alloy heat treated to 200,000 psi. From 
MIL-HDBK-5A, Feb 8 , 1966 : 

£ 

Elastic modulus in tension and compression, E = 29 x 10 psi 
Modulus of rigidity, G = 11 x 10 psi 
Therefore 

Poisson's ratio, v = .32 
Density, cu = O .283 #/in^ 

0.002 Compressive yield stress, a = 198,000 psi 

cy 



Total projectile weight: 



Warhead with filler 


9.74 lbs 


Fuse 


I .69 lbs 


Motor case 


1.99 Its 


Propellant 


3-59 Its 


Forward bulkhead 


0.40 lbs 


Nozzle plate 


I .89 lbs 


Fins 


0.31 lbs 


Fin Release 


0.10 lbs 


Obturator 


0.03 lbs 


Nozzle plug 


0.06 lbs 




19.80 lbs 



J 



forward = 11.43 lbs 



motor = 5.58 lbs 



aft =2.79 Its 



Length, L = 9-5 in 
Radius, r = 1.5 in 
Thickness, h = 0.08 in 
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DESCRIPTION OF THE EXTERNAL LOAD 



The external load on the projectile consists of two parts; l) the 

breech pressure at the base of the projectile, and 2) a lateral pressure 

on the sides of the case due to the breech pressure seeping beyond the 

obturator. A typical plot of the predicted barrel pressure as a function 

of time is shown in Fig. 1. According to NWC (Data sheet for 5 inch 

projectile, L. R. , 10-20-71) the base pressure on the projectile is 

0.9424 x breech pressure. Thus, a , the predicted maximum rigid body 

max 

acceleration of the projectile, is 

maximum „ 

force (26.9 ksi ) (0. 9424 )(tt) (1.5 in.) o QOf) 
max total 19.8 lbs/g ® 

mass 

where g is the acceleration due to gravity. 

The lateral pressure in the motor case occurs when the obturator 
cracks . Photographs of the firings show smoke coming from the gun 
barrel before the projectile emerges, thus establishing the existence 
of this pressure. No estimates of the magnitude of this pressure are 
available. Future NWC obturator designs will attempt to eliminate this 
pressure . 



DETERMINATION OF THE INTERNAL AXIAL FORCE 

In order to obtain an estimate of the transient internal axial force 

in the motor case, the projectile subjected to the transient base pressure 

shown in Figure 1 is idealized as a one dimensional steel structure with 

a force P sin 0 £ t £ .012, at the aft end. The forward end of 

o .012’ ’ 

the projectile is taken as either free or fixed, these two conditions 
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being the extreme possibilities. The belief was that this analysis 
would provide a bounded estimate of the transient axial force. The 
details of the analysis and the results are presented in Appendix A. 

The essential result of the analysis is that the rise time of the 
pressure is very slow in relation to the wave speed of the case, and 
hence the internal axial force at any location along the case at any 
time is equal to the mass forward of that location multiplied by the 
rigid body acceleration of the projectile at that time. Thus, there 
is neither an amplification nor a reduction in the internal axial 
force due to dynamic effects. 

Since the axial force in the case is dependent upon the distribution 
of the mass of the projectile, the question of the load path of the 
motor propellant becomes important. The propellant body force could 
either be transferred to a longitudinal shear load on the case or a 
compression load on the forward bulkhead of the aft end, depending 
to some extent upon the bonding condition between the propellant and 
the case. Since there was considerable uncertainty as to which path 
was the correct one, a two dimensional stress analysis was performed 
on the propellant for both the condition of bonding and of no bonding. 
The Rohm and Hass computer program AMG032A for the stress analysis of 
linear elastic, isotropic, axisymmetrically loaded, axisymmetric 
propellant gain with (or without) an orthotropic elastic motor case 
was used (Ref. l). The details and results are presented in Appendix 
B. The essential conclusion drawn from the results for the bonded 
propellant is that the maximum shear stress on the case (and in the 
bond) is approximately 300 psi and approximately 80 % of the total 
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propellant body force is carried by longitudinal shear loading on the 
case and 20$, by compression on the bulkhead at the aft end of the 
case. Of course, the unbonded propellant body force is carried 
entirely by a compression on the bulkhead. 

Taking the results of these analyses into consideration, the forces 
at the forward and aft ends of the motor case are 

^forward = ^ < 8 ’9°° ®> = ™ *ips 

P aft = 102 kips + ( 1>99 + g 8 X 3 ' 39 lbs )(8,900 g) 

= 145 kips 

for the bonded propellant (80$, carried by the case) and 

P aft = 102 kips + (^p)( 8 ’900 g) =120 kips 
for the unbonded propellant. The maximum stress due to these loads is 

°max = ( li|5 kips )/[(2 tt) (1.5 in.)(0.08 in.)] 

= 192 ksi 

which is just below a 

ey 

DETERMINATION OF THE DYNAMIC RESPONSE OF THE CASE TO THE TRANSIENT AXIAL 

FORCE 

Because the axial force is on the case for only a short period of 
time, the possibility exists that the case may not suffer large deforma- 
tions for loads that are larger than the critical load for static 
buckling, i.e., the loading maybe impulsive. This possibility is 
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investigated in Appendix C , where the critical static buckling mode of 

nt 

the case is subjected to the force P q sin - and the elastic 

transient response is computed. Two estimates of the natural frequency 

of the critical buckling mode are used, one without the propellant 

mass attached and one with the propellant mass attached. The essential 

result from this analysis is that if the maximum value of the applied 

force, P q , is equal to the critical load for static buckling the 

maximum magnitude of the buckled deflection is approximately 100-200 

times larger than the initial imperfection in the shape of the 

critical buckling mode. For the manufacturing tolerances in this 

program this is a large deflection and represents a severe buckling 

condition. Thus, the loading is not impulsive and the maximum axial 

force must be less than the static buckling load to prevent severe 
* 

buckling . 

DETERMINATION OF THE STATIC BUCKLING LOADS 
OF THE CASE 

The results of the preceding analyses indicate that the motor 
case can be designed using standard procedures for estimating the 
static buckling strength of shells. Consequently, the following static 
loading conditions were considered: 

(la) Uniform axial compression, only, P = Pf orward = p a ^ 

(lb) Linearly varying axial compression only, P = Pf orwar ^ + tt( 3 in)(x)q g 

* , 

A value of P equal to 90% of the critical load for static buckling 
causes an ampiification of approximately 10 which may also represent 
severe buckling. This depends upon the magnitude of the imperfection 
that exists in the buckling mode shape prior to firing. 
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(2) Uniform external radial pressure only, q r 

(3) Combined uniform axial compression and uniform radial 
pressure, P - P foMard = P aft , ^ 

where a and q are the radial and longitudinal uniform external 
^ s 

pressure components respectively, and x is the longitudinal coordinate 
(x=0 at the forward end) . 

The elastic buckling load for loading condition (la) was determined 

using the formula presented in Ref. 2, page 528, and the computer program for 

the geometrically nonlinear analysis of arbitrarily loaded shells 

of revolution developed by Ball (Ref. 3). The results are shown 

in Fig. 2, where cr cr is the critical axial buckling stress. The 

maximum predicted axial stress in the case due to the firing is 192 

ksi and is also indicated in Fig. 2. Since o cr is much larger than 

ct , the buckling is inelastic and the predicted elastic buckling 
cy 

loads are invalid. The inelastic buckling load for loading condition 
(la) was estimated using the charts from Ref. 2, page 704, and is also given 

-X- 

in Fig. 2. Note that this stress is less than the maximum predicted 
axial stress in the shell. As a consequence of this, the case will 
buckle inelastically due to the axial load alone . 

The elastic buckling load for loading condition (lb) was determined 
using the computer program of Ref. 3> and is shown in Fig. 2. For 
this condition all the propellant was assumed to be carried by a 
uniformly distributed shear stress on the motor case. The ratio of 



The result shown in Fig. 2 is a conservative estimate since NWC experimental 
results for the static buckling of a O.065 in. thick cylinder indicated 
a higher buckling stress. 
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axial stress at the forward end to the stress at the aft end was 



^ = *674, the same ratio as the total mass ahead of the 
forward end to the total mass ahead of the aft end of the case. The 
critical stress shown in Fig. 2 is the stress at the aft end. 

Buckling due to loading condition (2) is elastic. The results 
from Ref. 2, page 540, and the computer program (Ref. 3) are presented 
in Fig. 2, where is the critical external radial pressure. 

For elastic buckling due to loading condition (3), Ref. 2, page 
544, recommends 



/- + -SL- 1 

P q 
cr u cr 



The result from the computer program is given in Fig. 2. Note that 
this result agrees with the recommendation of Ref. 2. No results 
for the inelastic effects for loading condition (3) are readily 
available . A computer program that can treat inelastic buckling of 
cylinders under general loading must be used to obtain the inelastic 
buckling load under that loading condition. Hence, the stability of 
the shell along the dashed line is unknown at the present time. 

Also shown in Fig. 2 is the maximum axial stress possible i.e. 

200,000 ksi. Hence, the cross-hatching in Fig. 2 indicates the 
limiting values of the allowable loading. From this it can be seen 
that if the lateral pressure, due to leakage beyond the obturator, is 
larger than 2800 psi the case will buckle. 

The critical elastic buckling mode predicted by the computer 
program for all four loading conditions was essentially sin cos 29 
where 0 is the circumferential coordinate. 
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The details of the analyses using the formulas and charts in 
Ref. 2 are given in Appendix D. The computer results have been given 
to NWC personnel, Code 4573. 

An examination of Fig. 2 reveals that there are significant 
differences in the predicted elastic buckling load obtained from Ref. 2 
and from the computer program for loading conditions (la) and (2). 
Apparently, there is need for further study here. The shell considered 
here is not exactly a thin shell, and the formulas in Ref. 2 may not 
be valid for this shell. 



DESIGN TECHNIQUE 



The following plan is recommended as a general design technique : 

1. Determine the loads P f orwar( j 5 ^aft’ ^r 5 an< ^ ^ rom 4ke estimated 
weights and acceleration data and a propellant stress analysis. 

2. Select the thickness of the case such that the maximum stress is 
less than 198 ksi. 



3. Check to see if P a ^ and q g need to be changed. 

4. Determine if the shell is stable ( inelastic ially) under the loads 

^forward’ ^aft’ ^r 5 ^s ’ Note that q^ should include both the effect 
of the slumping propellant and the effect of the external pressure. 
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CONCLUSIONS 



Several analyses have been performed on the 3-inch motor case to 
determine the cause of the buckling failure. The results of these 
analyses are: 

1. The dynamic effects are insignificant and the case responds to the 
breech pressure as if the pressure were applied in a static 
sense, i.e. the maximum internal axial load in the case can 

be computed on the basis of the mass distribution and the 
maximum rigid body acceleration of the projectile. 

2. The loading is not impulsive, i.e. it is on the case for a 
sufficient length of time such that the shell will collapse if 
the applied load equals the static buckling load of the case. 

3. Approximately 80$ of the propellant inertial load is carried 
by shear along the case when the propellant is bonded. 

4. The case will buckle inelastically due to the axial load caused 
by the launch. 

5. If the obturator cracks, 10-20$ of the breech pressure acting 
as a lateral pressure is sufficient by itself to cause elastic 
buckling. 

6. The predicted buckling mode for axial load and for lateral 
pressure is the same and is essentially identical to the 
observed buckled shape of the launched projectile. 

7. The static buckling loads computed using the formulas in 
Ref. 2 are considerably lower than those predicted by the 
computer program of Ref. 3 • 



16 



REFERENCES 



1. Becker, Eric B. and Brisbane, John J., "Application of the Finite 

Element Method to Stress Analysis of Solid Propellant Rocket 
Grains, Vol. I and II," Rohm and Hass Company, Redstone 
Arsenal Research Division, Huntsville, Alabama, Report No. 5-76, 
Nov. 1965 and Jan. 1966. 

2. Baker, E. H., et al., "Shell Analysis Manual," MSA Contractor 

Report CR-912, April 1968. 

3. Ball, R. E., "A Geometrically Nonlinear Analysis of Arbitrarily 

Loaded Shells of Revolution," NASA Contractor Report NASA 
CR-909, Jan. 1968. 



17 



APPENDIX A 



Transient Axial Response 
Statement of the Problem 

In this section, the case is idealized as a one dimensional, 
one material bar with an input force function P(t) = P q sin ~ , 
where a is the duration of the force and P q the maximum intensity 
at t = tt/2. The time length of the input force was taken as .012 sec. 
The boundary condition at the aft end of the bar (x=0) was taken 
as free, that is the axial force in the bar is equal to the forcing 
function. Two boundary conditions were considered at the forward 
end of the bar (x=l); 

i) a free end with axial force equal to zero, and 

ii) a fixed end with displacement equal to zero. 

The actual boundary condition at the forward end is somewhere between 
these two idealized conditions. Hence the actual solution lies 
between the bounds obtained from the free-free and free-fixed problems. 

The physical, parameters taken in this problem, were the length 
of the bar i = 10 in . , and the Young ' s modulus for a steel bar 

£ O 

E = 30 x 10 lbs /in. . 




A-l 



Governing equation for axial response 




.2 d 2 u 




for moving bar 
0 for stationary bar 



2 






where u = displacement along x axis, x axis attached to the bar 
t = time 

.fp 

c = speed of wave propagation = | — 

E = Young's modulus of elasticity 
p = mass density 
Z = length 
A = area 

The initial conditions are 



u (x,0) = 0 
u t (x,0) = 0 



where u, denotes velocity • 

t or 



The boundary conditions are 
i) for free-free boundaries 



u x (0,t ) = G(t) 
u x (Z, t) = 0 



ii) for free-fixed boundaries 



u x (0,t) = G(t) 
u (A,t) = 0 



In equations (5) through (9) 



u 



X 



denotes — 
dx 




where A is the area of the bar. 



I. Axial response with free-free boundary conditions with 

P(t) = P sin — (10) 

o a ' ' 

Then 

P 0 . . 

G(t) = =£ sin — = P sin — (ll) 

EA a a 

We seek solution of equation (l), with initial conditions (3) and (4), 
and boundary conditions (5) and (6), with forcing function given by 
( 10 ). 

Solution by superposition. Consider solution u(x,t) in the form 

u(x,t) = v(x,t) + U(x,t) (12) 



Then (l) becomes 



tt 



+ U 



tt 



- c 



(v 



XX 



U ) 

XX 7 



P(t) 2 P . rrt 
= -c -r sm — 
p H A JL a 



^ (13) 



Letting 



f(x,t) = - U tt + c U 



2 „ _ P(t) 

xx p l A 



equation (13) becomes 



v, , - c 



V = f(x,t) 



tt XX 

Boundary conditions (5) and (6) become 

v (0,t) = G(t) - U v (0,t) 
v x u,t) = -U x (JL, t) 

and initial conditions (3) and (4) become 

v (x,0) = -U (x,0) 



v t (x,0) = -U t (x,0) 



(14) 

(15) 

(16) 

(17) 

(18) 



A- 3 



(19) 



U(x,t) is chosen so that equations (l6) and (17) become homogeneous 
boundary conditions in v. Here we take 

2 

U(x,t) = (x - |j) G(t) 

With equation (20) we find 

i) boundary conditions (l6) and (17) become 
v (0,t) = 0 

J\. . 0 

v x (X,t) = 0 

ii) initial conditions (l8) and (19) become, 

2 

v(x,0) = - U(x,0) = (x - ^-) G(0) 

But G(0) = 0, hence v(x,0) = 0 

2 

v t (x,0)= - U t (x,0) = (x - |j) G t (0) = Y(x) 

Since G (0) = P J 7(x) = ^ (x - |j) 
t 

v(x,0) = 0 



where 



v t (x,0) = 7(x) 



Y(x) - ^ (x - fj) 



( 20 ) 



( 21 ) 

( 22 ) 



( 23 ) 

( 24 ) 

( 25 ) 



The solution to the boundary value problem in v(x,t) defined by equations 
(15), (21) - (24) is sought . A solution in the form of a Fourier 
series is assumed, say 



v(x,t) = ^ v n ( t ) cos “ 



TTX 



(26) 



n=0= 



which satisfies boundary conditions (21) and (22). Functions f(x,t) 
and Y(x) are also expanded into their Fourier series. 



(27) 



where 



f(x,t) = ^ ** n (t) cos 



n=0 



f o (t) = i J d 5 



0 

a 



f n (t)= l J cos n > 0 



0 



(28) 



(29) 



and 



7(x) = £ 



... n tt x 

Y COS ;; — 

n A 



where 



n=0 

A 



’o'H ' f (?) d 5 



(30) 



(31) 



\ = § J 7 (?) cos d | n > 0 



0 



(32) 



Substituting expressions (26), (27) and (30) into equations (15)? 
(23), and (24) yields 

Z „ v r ^ v (t) .2 .. 

I “• W { ~ V- * ( ¥) *„<*> - f n<*>} ■ 0 



n=0 



dt 



I COS 2J1£ v n (0) - 0 



n=0 



£ i v (0) v 

I f-fr- - - 



0 



(33) 



(34) 



(35) 



n=0 



A- 5 



Equations (33) - (35) are satisfied if for each n, (n=0,l,2, . . . ) 



d 2 v n (t) 2 




v (t) 
n 



= f (t) 
n 



v (0) = 0 



dv (0) 
n 

dt 




( 36 ) 

( 37 ) 

( 38 ) 



Equations ( 36 ) - ( 38 ) may be numerically solved on the computer (using 
Runga-Kutta or some other scheme) for v n (t), for n=0,l,2,..., where 
the tilda denotes the numerical approximation solution. Hence for a 
convergent series, 

00 00 

v(x,t) = ^ v n (t)cos v n (t)cos 2 -y^=v(x,t) ( 39 ) 

n=0 n=0 

Once v(x,t) is obtained, we return to equation (12) to find 

u(x,t) = v(x,t) + U(x,t) =- v(x,t) + U(x,t) = u(x,t) (40) 

Hence , 

N 

u(x,t) =“ ^ v n (t) cos + U(x,t ) (4l) 

n=0 

The force is given by 

F = “ £ - {- l ¥ ^ + I (x,t) } » <* 2 > 

n=l 

Note that F does not depend upon the n=0, or rigid body term. 
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Numerical solution of 



d 2 v (t) ? 2 



v (0) = 0 
n 



f (t) 

n 



n = 1,2, 



dv (0) 
n 

dt 



= Y_ 



n 



( 36 ) 



( 37 ) 



( 38 ) 



Recall eq. (29) for n > 0 

l 



f n (t) = | J f(5,t) cos d ? = !/ [(-? + i) G tt - niS d| 



= I k + f- , \g cos 5 sin - 

A L tt l /nnY~ A ^nrr^ * 

. nn? , § . n tt 5 \ 

m — 7-* + — T 2 — r- sm — 7^- f 
i _ - /nir\ l J 



¥ + cos 2-4 

1 2i n ^ 2 



<f) £ 



«(f) 3 



f) 



-(-« 



+ 44 A sin SiS f 



r u 7S y hs “- 1 " 1 



Therefore 



or 



f (t) 
n 



M c 

2 2 u tt 

n rr 



21 ^ . TTt 

~2 2 Psln — 
a n 



xr . TT *fc 

f = K sm 

n n a 



where K = 
n 



2.1 P 

2 2 
n a 



( 44 ) 

(45) 

( 46 ) 
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Let P n = j— , then eqs. (36) to (38) become, 



n 



d P (t) « . N 2 . 

" + c 2 (H!) e = sin 25. 

dt 2 Vi/ n or 



Mo) = 0 

n v ' 



dp Y 

__E ( 0 ) = -2. 
dt v ; K 
n 



Recall equations (32) and (25), 

Z 

\ - ! J ' f (5) =0= ^ d ? 



0 



. .2 f" Si (, . £) cos 2-Si . 2nP ' -i 2 \ 
n X J a ^ 2r COS jl sM 2 2/ 

n tt 



Y 

= + - = + 3-82 x 10" 3 
K n 11 

for a = .012 sec. 

Hence we solve for P n (t) from the following equations 

.2 



11 + 0 2 (H!) e = sin Si 
, . 2 \ Z 1 K n a 



ij 
dt‘ 

P (0) = 0 
n v 



dp (0) 

— — = - = 3.82 x 10 _:5 for a = .012 
at tt 

After obtaining P (t) from the computer we find 



n 



v (t) = K p = - —rT~v 
n n n 2 2 

n a 



Z P P (t) 
n v 



(*+7) 

m 

( 49 ) 



(50) 

(51) 



(52) 

(53) 

(5^) 
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(55) 



Finally then the force per unit P is 



N 



vL(t) 



F _ / V nn n v ' . n ttx , 1 3U(x,t) \ 

f rl t T sin — r + 5 j “ 

n-l 



Therefore 



N 



F V 2 tt a nrrx , x % . rrt 

P“ = L ~ P n (t) sin ~r + (1 “ sln ~ 



(56) 



n=l 



a n 



The latter term is the internal force due to the rigid body acceleration. 
The terms voider the summation sign introduce the transient dynamics 
due to the elasticity of the bar. 



II. Axial response with free-fixed boundary conditions with 

P(t) = P o sin 752 ( 10 > 

A solution of equation (l), with initial conditions (3) and (4), and 
boundary conditions (7) and (8), with forcing function (10) is sought. 
The method proceeds as before, except here we take 



U(x,t) = (x-JL) G(t) (57) 

Boundary conditions on v(x,t) become 

v x (0,t) = 0 (58) 

V u,t) = 0 (59) 

The governing equation in v(x,t) is 

v tt " ° 2 v xx = ( 6o ) 



A-9 



where 



f(x,t) - - u tt + c 2 Uxx ( 61 ) 

= - (x-1) G tt 

The initial conditions are 

v(x,0) = 0 ( 62 ) 

v t (x,0) = (x-X) G t (0) = Y(x) (63) 

where 

Y(x) = (x-X) (64) 

ot 

Assume the Fourier expansion 

00 

v(x,t) = £ v n (t) cos ( 65 ) 

n=l,3, • • • 



which satisfies boundary conditions ( 58 ) and 
and 7(x) are expanded into Fourier series, 

00 



f(x,t) 



I f n (t) cos 
n=l , 3 > • • • 



(59). 



As before f(x,t) 



( 66 ) 



21 

V*) ■ h, J f(5 > t) « 

0 

2SL 

f n (t) = J J f(?,t) cos d| n > 0 
0 

00 

tU) = £ Y n cos Sjp 

n=l,3, . . • 

21 

\ = j J Y(S) cos d? n > 0 

0 

21 

% - h J ' f (s) cos ^ a 5 
0 



( 67 ) 

( 68 ) 

( 69 ) 

(70) 
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Substituting eqs. (65) - (70) into eq. (60), (62) and (63) yields 



d 2 v (t) 



I cos { TJ - + ° 2 fl* V n (t) - f n (t) } ‘ 0 

n=l 5 3 5 • • • 



dt 



1 



cos 



n TTX 

21 w n 



v (0) = 0 



n=l,3, • 



“ r dv 0 , 

U 2 i 1 dt n J 



0 



n=l,3, • • • 



Again we seek the solution to the equations 
2 

d v „ , s 2 



n , 2 / rm \ „ 

l? + ° V n =f n 



( 71 ) 



(72) 



(73) 



(74) 



v (0) = 0 
n' ' 



dv (0) 
n 



, = v 

dt n 



(75) 



(76) 



for n = 1,3,5,... Numerical solution for a finite number of n = 1,3,5, 



N yields the approximate solution, v (t), hence 



n 



N 



v(x,t) = ^ v n (t) cos 

n=l,3,... 



x 



Then 



N 



u(x,t ) =- ^ v n (t) cos + (x-4) G(t) 

n=l,3j • • • 



Finally 



Su(x,t) 

Sx 



N 



_ V nTT ~ / . \ . n TT X _ / , \ 

■ l 21 v n (t) sin ~Z T + G < t > 



n=l,3> • 



(77) 



(78) 



(79) 
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Numerical solution of 



d 2 v 



dt' 



n , 2 

? + c 



rm V 
21 ) n 



v = f 



n 



v (0) = 0 
n v ' 



dv (0) 
_n v ' 

dt 



= Y„ 



n 



for n = 1,3,5: 



Calculation of f (t) and Y for n > 0 
n v n 

21 



21 



f n (t) = T J cos ^ « - T I d-S> G tt cos ^ (80) 



— f (t) = 



G tt n 



0 

8i 



2 2 
n tt 



for n = 2,4,6, . . 
for n =1,3,5,.. 



(81) 



\ = I J *<?> COS ^ d5 . 1 J I !2_ (S-i) cos ^ d| (82) 



0 



V = nP 

T n .012 



0 for n even 
- -8 & 



2 2 
n tt 



for n odd 



(83) 



Hence we seek the solution of 



d 2 v (t) 



n 



+ c 



dt 






V = 

n 



-8 X P 



2 2 
a n 



sm 



TT t 

.012 



v (0) = 0 
n 



dv 

at 2, (0) = 2 2 
n tt 



-8i 



rr P 
a 



for n=l,3,5,. 



(84) 

(85) 

( 86 ) 



Let K = Jr - o j and define 6 = v /K , 

n 2 2’ n n' n 



2 2 
a n 



(87) 
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then the above equations become, 




n 



2 

+ c 




sin 



TT t 

a 



(88) 



p (o) = o 
n x / 



n = 1,3,5, • . • 



(89) 



dp _o 

-r- 7 ^- (0) = — = + 3*82 x lO - ' 5 for a = .012 sec 

dt ' TT 



(90) 



APter B is determined, we have 
n 



v (t) = K P = ~ 8 9 P (t) 

n n n (.012) 2 n 2 n 



n=l ,3,5. 



(91) 



and 



F_ 

PO 



N 

= V 

L-> 

n=l , 3 ? • • • 



— P (t) sin 
a n 



n rrx 
21 



+ sin 



Tft 

a 



(92) 



for the solution of a free-fixed bar. The latter term represents the 
static response to the end force. The terms under the summation 
sign introduce the dynamic effects. 



III. Results 

The time history F /P at x/ i = 0,0.25,0.5,0.75 and 1 are given 
in Figs. A1-A5 for the free-free bar with a = .012 sec. Figs. A6-A10 
show the corresponding values of F/P q for the free-fixed bar. Figure 
All shows the history of f/P q at x/JL = 0.5 for a = 0.0012 and Fig. A12 
shows the history at x/JL = 0.5 for a = 0.012 and F = 3 x 10 psi for 
the free-free bar. These figures are from a plotter. The solution 
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t F /p. 



•H 

a 
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Fia A 2. 




a-i6 
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FREE - FR.6& c**.0ia 



was plotted after each new term in the series was added. The solution 

with the contribution from U was plotted last. Hence the small 

wiggles indicate the dynamic contributions. Figure A13 shows F/P 

o 

as a function of x/X for several values of t for a = 0.012 sec. for 
both bars. 

IV. Summary of Results 

The solutions of the boundary value problem show that for a forcing 
function of time duration .012 seconds, the force appears as a standing 
wave. That is, the time length of the force is so much larger than 
the time of wave propagation that the dynamic effect is negligible. 
Decreasing the time of the force pulse to .0012 seconds showed a 
noticeable increase in the dynamic effect. However even at this 
shorter time pulse, the internal forces are but slightly larger than 
those for a = .012 seconds. A further decrease of a to .00012 seconds 
show that the dynamic forces are the same order as the static result 
and there is an amplification of the internal force in the bar. 

With the time pulse at a = .012 seconds, Young's modulus of the bar 
was reduced from 30 x 10^ psi to 3 x 10^ psi to determine the effect 
of "softening" the bar material due to inelastic effects. It can be 
observed in Fig. A12 that although a noticeable increase in the dynamic 
effect occurs the internal force is not significantly altered. 

The conclusion that the dynamic effects are insignificant can be 
supported by an examination of equation (52). For each value of n, 
equation (52) is equivalent to a single degree of freedom system with 
a harmonic forcing function of frequency rr/a rad/sec or 260 rad/sec 
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for a - 0.012 sec. The natural frequency of the single degree of 
freedom system is cnrr /i. For l = 10 and c = 2.02 x 10^ in/sec, the 
lowest natural frequency is 63,500 rad/sec. Hence, the ratio of the 
forcing frequency to the lowest natural frequency is 0.004l. The 
magnification factor for this frequency ratio is essentially one, 
i.e. the case response in a static sense. 
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APPENDIX B 



Two-Dimensional Stress Analysis of the Motor Propellant 

The results for the propellant -case interface stresses from the 
Rohm and Hass computer program AMI032A (Ref. l) are given in Fig. B-l 
for E = 5 x 10 psi and a body force of 480 g, which is equivalent to 
an acceleration of 8000 g of the propellant grain whose density is 

Q ■X' 

0.06 lb/in . Also shown in Fig. B-l is the location and the magnitude 
of the maximum principal tensile stress. In this analysis the motor 
case was assumed to be rigid. The nodes of the finite elements used 
to model the propellant are denoted in Fig. B-l by I, J, where I 
defines the radial location and J the axial location as shown in Fig. 

B-l. The stresses predicted by the program are at the center of each 
element. Hence, the stresses shown in Fig. B-l are actually located 
0.125 in. from the propellant-case interface. A plot of the shear 
stress from the inner face of the propellant to the outer face is 
shown in Fig. B-2 for J = 14. Note that the outer face propellant 
shear stress is approximately 10$ larger than the stress given in 
Fig. B-l at J = 14. 

The total vertical force provided by the longitudinal shear stress 
along the case and the axial compression stress on the bulkhead at 
the aft end is 

•x- 

The computer program output consists of the coordinates and displacements 
of each nodal point and the stresses and strains for each continuum 
element . 



B-l 



■ &oc> -^/po -goo 



1 = / 




0"l Prop* I \an~h 



arh Tnj~er-fac.c:^ Boodec 
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Reaction = (0 ‘ 25 in *)( 2 57 + 273 + ••• 58 + 19 psi) (1 .1 ) (tt) ( 3 in.) 

+ (0.25 in . ) [ ( 930 psi) (1.25 in.) + (787 psi)(l.75 in.) 

+ (832 psi) (2 .25 in.) + (731 psi) (2.75 in.)psi](TT) 

= 22.9 kips (shear) + 5-0 kips (compression) = 27.9 kips 

when the shear stresses shown in Fig. B-l are increased by 10$ to 
extrapolate to the interface stresses . The total force applied to 
the propellant is 

P applied = (0,06 lb /in. 3 )(1.5 2 - 0.5 2 )(in. 2 )(9-5 in. ) ( tt) ( 8,000 g) 

= 28.6 kips 

Thus, the stress analysis appears to provide total equilibrium. The 

fact that P , . is slightly less than P . . , is due to the fact 
reaction applied 

that the compression stresses at the propellant -bulkhead interface 

are somewhat larger than the values shown in Fig. B-l, which are 

at the center of the bottom elements 0.125 in. above the bulkhead. 

From these results 0 *q = 82$ of the total propellant inertial 

c. ( .y 

load is carried by longitudinal shear stress on the motor case. 

4 

The same problem was analyzed with E = 5 x 10 psi. The results 
were essentially identical. 

5 

The stresses for the unbonded propellant with E = 5 x 10 psi are 
shown in Fig. B-3. The total vertical force is carried entirely by 
compression on the bulkhead at the aft end and is equal to 

Preaction = ( °* 25 in )£( 4 > 418 psi)(l.25 in.) + ( 4 , 4 90 psi)(l.75 in.) 
+ (4,516 psi)(2.25 in.) + (4,526 psi)(2.75 in.)](rr) 

= 28.2 kips 

The computer output for these three runs have been given to NWC 
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APPENDIX C 



Transient lateral Motion 
Statement of the Problem 

The object of this investigation was to determine the transient 
lateral behavior of an axially loaded bar with initial deformation 
w 1 (x). The bar is idealized as simply supported at both ends, and 
is subjected to an axial compression P(t) = P q sin where 0 / = .012 
seconds . The behavior of the bar was considered for various ratios 
of P q , the maximum intensity of thrust, to P g , the Euler buckling 
load for a simply supported column. 



C-l 



w 4 

L 



P(t) 






B .C . ' s simply supported 



' vf 

4 y ; w 

w 1 - initial imperfection, assume w 1 



- y 



i . n TTX 
w sm — t— 
n i 



( 1 ) 



n=l 



- bending deformation, assume 



- I 



n TTX 



f n(t) sin t 



( 2 ) 



n=l 




) M. = 0: M - (M + dM) + (m + (V + dV) dx + P dw = 0 

u A at 



™ + v + ^=o 

dx dx 



^ V = 0: V - (V + dV) - m ^-| dx = 0 — ^ = m 

at 



-2 


(3) 


a w 
at 2 


(4) 
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Recall the moment curvature relation, 



M(x) = - El w 



XX 



Using (3) and (5), equation (4) becomes 

.4 d *2 -.2 

,,, d w _ d w S w 
-El r- - P — =■ = m — 

Sx Sx St 



With (l), (2) and ( 7 ), equation (6) becomes 



( 5 ) 



(6) 



- ^ El (££) 4 f n (t) sin + P(t) 

n=l 



I 

n=l 



(— ) 
k i ’ 



(w 1 
v n 



f (t)) 
n v ' ' 



sm 



n TTX 

A 



= m I sin ( 8 ) 

n=l 



which may be written 

I [ EI (f ^ <T > 2 <"n 

n=l 

n = 1,2, . . . 

Equation (9) is satisfied if each of 



is zero, i.e. 








EI ( T )4 f n (t) 


- p <*> <T> 2 K + f „ 


(t)) + m f^(t) = 0 


(10) 


with I.C.'s 


w(x,0) = w i (x) 




(11) 




|jr (x,0) = v(x) 




(12) 



+ i* n (t)) + m f^(t)J sin — ji = 0 

(9) 

the coefficients of sin n TT . X 
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Using 



P = 
e 



2 2 

n n El 



2 2 
n tt 

CO — 7Z 

n „2 




equation (10) becomes 



f n (t) * ® n a 



*»<*> 



P(t) 2 i 

-1 co w 
P n n 

e 

n = 1,2, 



with initial conditions, 



(13) 

(14) 



(15) 





f (0) = w 1 
n v ' n 


( 16 ) 




f (0) = 0 
n v ' 


(17) 


where 


i 






i P i/ \ . nrrx, 

w n =j w (x) sin — -jj— dx 


( 18 ) 



o 



For an impulse load in the form, 



P(t) ■ P o Sin T55T 



(19) 



and introducing the amplification factor 

f (t) 

^(t) = -V- 

w 

n 



( 20 ) 



the previous system of equations become 5 

g^t) + <4 [l - j? sin -g|] = j? 

e e 

g^o) = 1 



ut 2 
sin — - ■ co 
.012 n 



g n (0) = 0 
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( 21 ) 



This equation is also valid for a cylindrical shell subjected to 

an axial load or an external pressure where u) n is the natural frequency 

til til 

in the n mode and P g is the buckling load in the n mode. 

The previous system was solved numerically on the computer using 

a Runga-Kutta integration method to yield the results shown in Figs. 

Cl-CQ. Several values of P /P were considered for both cu = 6,000 rad/sec 

o' e n ' 

and 10,000 rad/sec. 

Summary of Results 



Amplification Factor |g(t) 


P /P 
o' e 


0 ) = 6000 rad/sec 


cd = 10,000 rad/sec 


.75 


3 


4 


.90 


10 


10 


1.0 


100 


180 


1.1 


600 


220,000 


1.2 


450,000 


* 

CO 


1.3 


* 


* 


CO 


CO 



j. , 

00 defined as greater than 1 x 10 [maximum g(t) possible 
from computer integration algorithm.] 



Estimation of the Natural Frequency 

. -£ ll a a + %) 4 

n a K m k 2 a 2 



m 



= mass/unit area 



From "Tables of Frequencies and ivfodes of Free Vibration of Infinitely 
Long Cylindrical Shells," by M. Barron and H. Bleich, J. Appl. Mech, 
Vol . 21, 1954 , pp. 178-184. 
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where 



a = 1.5 in. 

G = 11 x 10^ psi 

m = 7.43 x i 0 ' 1 * ^ 4 ° x .08 in. = 5-94 x 10 -5 ^~ s !! c 

in in^ 

when the mass and elasticity of the propellant is neglected. 

For w = sin sin 20* (— = 6 ) 

AJ 3 * 

K = 0.09178 
and hence 

cu = 9.6 x 

Adding the mass of the propellant gives 

m « 5.94 x 10' 5 + I .38 x 10’ U x tt x 1.5 2 /(3 tt) = 1.62 x 10 _l+ ■ ^" S( ! c 

in J 



v2 

? = 2331 



10 ^ rad/sec 



Hence 



cu = 5.8 x 10 



3 



rad/sec 
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APPENDIX D 



Computation of the Static Buckling Loads 



Axial Compression Only (Ref. 2, page 528) 



cr _ r Eh 
7] c r 



where o cr is the design - allowable buckling stress, i.e., nine out of 



ten shells will not buckle at ct . For elastic buckling T] = 1. From 
Fig. 3-23-1, page 530 of Ref. 2. 



C = 0.23 
c 



for E “ 53? 3 l8 ' 8# and 



- l’W) 1 * 

rh 



( 9 .^) 2 ( 1 -. 32 2 ) 1//2 

(l.5)(0.08) 



= 715 



Therefore 



cr 

H 



= (0.23) (29xl0 3 ksi) (0.08 in.)/(l.5 in.) 



= 356 ksi 



which is larger than ct - Hence the buckling is inelastic, i. 

cy 



7 ] < 1 . 



This value of r/h is beyond the end of the curve, and hence the value 
selected for C c is an estimated one. This indicates that perhaps this 
value is invalid for this small r/h ratio. 
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Figure 3.62-8, page 704 of Ref. 2, gives the critical inelastic 
buckling stress o cr as a function of c cr /Tl for elastic buckling 
for alloy steel 4130, heat treated to 180 ksi. From curve E^, 

a = 162 ksi 
cr 

This stress is based upon a 0.002 compressive yield stress of 164.5 ksi, 
which is considerably below the stresses of 179 ksi and 198 ksi given 
in MTL-HDBK-5A and for 4130 alloy steel heat treated to 180 ksi and 
200 ksi respectively. Increasing the predicted critical inelastic 
buckling stress to account for the higher heat treatment leads to 

a cr « (162 ksi) (|||) = 179 ksi 

when the MIL-HDBK stresses are used. 

A static buckling test was conducted at NWC on an axially loaded 
3 in. cylinder of 4130 alloy steel, heat treated with h = O.065 in. 

For this cylinder, 

= 290 ksi 
71 

and hence from Fig. 3-62-8 of Ref. 2 the inelastic buckling stress is 

a = 159 ksi 
cr 

so that 

a cr ~ (159 ksi) (j||) = 176 ksi 

From the experimental results, the axial force at buckling was 121.3 
kips . Hence 

<°cr>exp - (ni^KX in.) * *"-5 *.i 
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which is considerably above the predicted value of 176 ksi and is 

essentially equal to a . This very high value is somewhat surprising. 

cy 

Nevertheless, in light of this the predicted inelastic buckling load 
of 179 ksi for the 0.08 in. case must be considered low (conservative). 
External Radial Pressure Only (Ref. 2, page 5*+0) 

cr 2 2 

CP _ y TT E /Il\ 

~ = *12(1-V 2 ) ^ 



From Fig. 3*23-7 of Ref. 2, lateral pressure only, 



K = 24 
P 

for Z = 715 * Therefore 



cr (tt) 2 (29xlo\si)(O.Q8 in.) 2 

T| (12) (1-.32 2 ) (9*5 in.) 



= 45 ksi 



which is below a . 

cy 

lateral pressure is 



Hence , 



a h 



cr 

r 



this buckling is elastic, and the critical 



(45 ksi) 



(0.08 in. ) 
(1*5 in.) 



= 2.40 ksi 
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